Experimental Consequences of the S-wave cos^k^) ■ cos{ky) Superconductivity in the 

Iron-Pnictides 
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The experimental consequences ol different order parameters in iron-based superconductors are 
theoretically analyzed. We consider both nodeless and nodal order parameters, with an emphasis 
on the cos(fci,) ■ cos(fcy) nodeless order parameter recently derived by two of We analyze the 
effect of this order parameter on the spectral function, density of states, tunneling differential con- 
ductance, penetration depth, and the NMR spin relaxation time. This extended s-wave symmetry 
has line-zeroes in between the electron and hole pockets, but they do not intersect the two Fermi 
surfaces for moderate doping, and the superconductor is fully gapped. However, this suggests sev- 
eral quantitative tests: the exponential decay of the penetration depth weakens and the density of 
states reveals a smaller gap upon electron or hole doping. Moreover, the cos(fca;) ■ cos(fcy) supercon- 
ducting gap is largest on the smallest (hole) Fermi surface. For the 1/Ti NMR spin relaxation rate, 
the inter-band contribution is consistent with the current experimental results, including a (non- 
universal) behavior and the absence of a coherence peak. However, the intra-band contribution 
is considerably larger than the inter-band contributions and still exhibits a small enhancement in 
the NMR spin relaxation rate right below Tc in the clean limit. 



Introduction - The recent discovery of iron-based su- 
perconductors with a transition temperature as high as 
55K has stimulated a flurry of experimental and theoret- 
ical activity^i'^i^i^iii^i^iiSiiii. However, a conclusive ob- 
servation of the pairing symmetry still remains elusive, 
with both nodal and nodeless order parameters reported 
in experimental observations. 

Numerical and analytic research suggests that the an- 
tiferromagnetic exchange coupling between Fe sites is 
strongi2,ii^iii. Owing to As-mediated hopping, antifer- 
romagnetic exchange exists not only between the near- 
est neighbor (NN) Fe sites, but also between next near- 
est neighbor (NNN) sites. Moreover, the NNN coupling 
strength J2 is stronger than the NN coupling strength 
Ji. The Ji — J2 model produces half-filled magnetic 
physics consistent with experimental neutron data—-. A 
nematic magnetic phase transition has been predicted in 
this modeli^iii, consistent with the experimental obser- 
vation of a structural transition preceding the spin den- 
sity wave (SDW) formation. This model suffers, however, 
from an important deficiency - it is an insulator whereas 
the real material is an, albeit bad, metal. We, however, 
believe that the spin-spin interaction insight is important 
to the physics of the iron-pnictides. 

In a recent paper—, two of us added electron itineracy 
to the problem and studied a t — Ji — J2 model with- 
out band renormalization. We found that the singlet- 
forming Ji — J2 interaction gives rise to four possible 
pairing symmetries: cos{kx) ± cos(fcj,), sin(fc2;) • sin(fcj,) 
and cos{kx) ■ cos(fey). The last two are strongly pre- 
ferred from an interaction standpoint when J2 > Ji, 
but only cos(fc2,) • cos(fcj,) matches the symmetry of 
the iron-pnictide Fermi surface: it is maximal around 
(0, 0), (tt, 0), (0, tt), (tt, tt) - the location of the Fermi sur- 
faces in the unfolded one-iron per site Brilloiun zone. 



Although we used a specific, two band model for our 
calculation^, our results are completely independent of 
any model, as long as the dominating interaction is next- 
nearest neighbor J2 and the Fermi surfaces are located 
close to the aforementioned spots in the Brillouin zone. 
Some order parameters (such as dxy = sin{kx) ■ sin^ky) 
and others) mismatch the Fermi surface symmetry and 
can be discarded. We note that cos{kx) ■ cos(fcy) changes 
sign between the electron and hole pockets in the Bril- 
louin zone. In this sense, it resembles the order param- 
eter proposed by Mazin through weak-coupling general 
arguments^^. At moderate doping, our gap is isotropic 
within the same Fermi surface, while changing sign be- 
tween electrons and hole pockets, but at relatively high 
doping cos{kx) ■ cos(fcy) exhibits some anisotropy even 
within the same Fermi surface. 

Neutron measurements have found antiferromagnetic 
stripe order of Fe moments ranging from in 
NaOFeAsiS, and O.Se/is^o in LaOFeAs to O.S^is in 
CeOFeAs2i and SrFe2As2^. A magnetic moment of 
0.8^B is fully consistent with a purely localized spin-one 
Heisenberg model. While a magnetic moment of O.^fiB is 
smaller than what is expected in a purely localized spin- 
one system, it is rather larger than what can be obtained 
in a truly weak coupling theory. We point out that, 
due to imperfect nesting, weak coupling theory requires 
large values of U/t ~ 4 to explain even small magnetic 
moments (< 0.2fiB), clearly outside the weak-coupling 
limit^^. Considering these facts, together with the rather 
high resistivity of the iron-pnictides, we find that the ex- 
perimental evidence paints a picture of the iron-pnictides 
as being at moderate interaction couplings. Thus, mod- 
erate to strong coupling models can provide an accurate 
qualitative description of the observed phenomena. In 
fact, the t — Ji — J2 model predicts the right physics of 
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the parent state SDW as well as the cos(A:2;)-cos(fcj,) order 
parameter. 

In this paper we focus on the experimental properties 
of several superconducting order parameters proposed 
in the iron-pnictides, with particular emphasis on the 
cos(fca;) •cos(A;y) order parameter. We look at a simphfied 
two-band superconducting model and obtain the spectral 
function, density of states, tunneling differential conduc- 
tance, penetration depth and NMR spin relaxation time. 
We stress the important point that the cos(fc2:) • cos(fcj,) 
order parameter features lines of zeroes at {±n/2, ky) 
and (fca;,±7r/2), as in Fig. [T] (obviously, irrespective of 
its harmonic form, any order parameter changing sign 
between the electron and hole Fermi surfaces must have 
zero-lines). Thus, at low doping, the hole and electron 
Fermi pockets are far away from the zero lines of the 
order parameter and the superconductivity is nodeless. 

Close to half-filling, we find that the cos{kx) ■ cos(fcj,) 
order parameter exhibits an exponentially decaying 
SX{T) = \{T) - A(0), where A(T) is the penetration 
depth at temperature T, as expected for a nodeless su- 
perconductor. However, upon doping, the gap on the 
Fermi surface varies in magnitude: for electron doping, 
the gap decreases on the electron pocket and increases 
on the hole pocket. The penetration depth is sensitive 
to the smallest gap in the system, and hence exhibits 
a weakened exponential decay upon doping. This could 
explain the conflicting values of the gap parameters ob- 
tained by fitting the penetration depth experiments to 
the BCS exponential form— i^^i^. In the unlikely event 
that the system remains superconducting at very large 
doping, then the Fermi surfaces will cross the line of ze- 
ros of cos(fc:E) • cos(fcy) at around 35% doping, and cause 
SX{T) to become linearly dependent on T. 

We also calculate the NMR spin relaxation rate 1/Ti 
of the bare superconductor and find that it factorizes 
into inter- and intra-band contributions. While, for the 
cos(fca;) • cos(fcj,) order parameter, the inter-band contri- 
bution to the NMR spin relaxation rate does not exhibit 
a coherence peak, the intra-band contribution is larger 
than the inter-band contribution and still exhibits an 
enhancement right below Tc owing to its fully gapped 
s-wave nature. Adding the two contributions we find 
that, although the coherence peak for cos{kx) ■ cos(fcj,) 
is smaller than for a sign-preserving gap such as, for ex- 
ample, |cos(fca;) • cos(fcj,)|, it is still present due to the 
intra-band contribution. The coherence peak can be 
strongly reduced if the intra-band scattering is stronger 
than inter-band scattering or if the samples are strongly 
disordered. If the As structure factor A{q) is taken into 
account, the inter-band contribution is severely reduced 
due to the fact that A{q) = cos{qx/2) cos(gy/2) is zero 
close to the wavevector difference between the electron 
and hole Fermi surfaces: q = (±7r, 0), (0, ±7r). The As 
structure factor also reduces the overall coherence peak 
by smearing the intra-band contribution. 

Model - We approximate the typical iron-based ma- 
terial by a two-dimensional square lattice of Fe atoms. 
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FIG. 1: Schematic diagram of the Fermi surfaces in the 
iron-pnictides at half-filling in the unfolded Brillouin zone 
— TT < kx < n, — TT < ky < n. The dashed (red) and solid 
(blue) curves correspond to the hole and electron Fermi sur- 
faces, respectively. The dashed lines mark the nodal lines 
at (±7r/2,fcy) and (fca;,±7r/2) for the cos(fca;) ■ cos(fcy) order 
parameter proposed in Ref. QJ. 
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FIG. 2; Behavior of the spectral function A(k, w) in an in- 
terval about the Fermi energy (—0.02 < ui < 0.02) across 
the unfolded Brillouin zone — tt < fc^, < tt, — tt < ky < n 
for gap parameter Ao = 0.1. Panels (a) and (b) depict the 
Sx2^y2 and dxy order parameters, respectively, both at chem- 
ical potential p = 1.6. The Sx-2y2 order parameter is shown in 
panels (c) and (d) for the higher electron doping values /j. — 2 
and jj, = 2.2, respectively. For these high doping values, the 
s^2y2 superconductor has become nodal. The lighter regions 
illustrate the ungapped portions of the Fermi surface. 
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since the superconductivity has been shown to be asso- 
ciated with the FeAs layer. To capture the degeneracy 
of the dxz and dyz orbitals on the Fe atorns, ^we use the 
two-orbital per site model proposed by Ref. 23. Although 
this description is only valid in the case of an unphysi- 
cally large crystal field splitting, we particularize to this 
model for analytic simplicity. The kinetic part of the 
Hamiltonian is written: 



with = — fJ-, = £y ~ fJ-, and we have used the four- 



ea;i/(k) 



ea;y(k) 



(1) 



Here, VL = i^i.: 



) is the creation operator for 



spin a electrons in the two orbitals (1,2) = {dxz,dyz), fJ- 
is the chemical potential, and the matrix elements are 

ea;(k) = ~2ti cos fcx ~ 2^2 cos ky — At^ cos cos ky 
ej,(k) = —2t2 cos kx — 2ti cos ky — 4^3 cos kx cos ky 

ea;y (k) = — 4^4 sin /cj; sin fcy (2) 

While Eq. ^ is only a simplified version of the true 
band structure of the material, it produces Fermi pock- 
ets that resemble those predicted by density functional 
theory (see Fig. [T|). The eigenvalues of ([T]) are: 



E± = e+ - fi ± 



xy 



(3) 



where t± — {tx^ty) /2. In the following, we take ti — —1, 
t2 — 1.3, and ^3 = t4 = —0.85. The undoped compound, 
where there are two electrons per site, corresponds to 
/I = 1.54. 

We now assume that the interacting part of the Hamil- 
tonian induces singlet pairing between electrons within 
each orbital, but we make no further assumptions about 
the form of the interaction or the pairing mechanism. 
Then we introduce pairing gaps A1.2 for each orbital 
and we write down the mean-field effective Hamiltonian 
i7(Ai, A2) = 2k*(k)^5(k)*(k), where 



B(k) 




Ai(k) 

-e.(k) 



— txy{k.) 



^xy{k) 

-ea;y(k) 

W A2(k) 

A^(k) -e,(k) 



(4) 



component spinor ^'(k) = (ci.k. 



C2,k,T,C2 k, 



...... V' 

We neglect inter-orbital pairing in order to make the 
problem analytically tractable. This is also reasonable 
because two of us proved in Ref. Uj that, at least for the 
case of the t — Ji~ J2 model (and hence for the most im- 
portant gap we will be focusing on - cos(fca;) •cos(fcj,)), the 
inter-orbital pairing expectation value is negligible even 
in the case of strong Hund's rule coupling. 

The symmetry of the superconducting order param- 
eter A(k) has two possible d-wave types^: dx2-y2 ~ 
Ao(cos kx — cos ky) and dxy ^ Aq sin kx sin ky, and three 
possible s-wave types^: Sx^+y^ ^ Ao(cosfc2: -I- cosfcy), 
Sx^y^ ~ Aq cos kx cos ky, as well as the constant gap (sq) 
which is not allowed in the t — Ji — J2 model but can 
obviously appear in other interacting models. The C4 
symmetry of the underlying lattice maps kx ^ ky and 
dxz dyz- Hence for all the pairing symmetries de- 
scribed above we have /S.i{kx,ky) = iS.2{ky,kx), except 
for (ia.2_j,2 where Ai (fc;j,, fcj,) = —l!^2{kx,ky)XT\ie dx2_y2, 
dxy and Sx^+y'i pairing symmetries are nodal while the 
other pairing symmetries are nodeless. We now proceed 
to analyzing the experimental consequences of these pair- 
ing symmetries. 

Spectral Function, Density of States and Tunneling 
Differential Conductance - The single-particle density 
of states (DOS) can be written as: 



AA(a;)^^^(k,c.) 

k 

= V3[^ll(k,c^ + ^<5) + g33(k,c^ + ^<5)] (5) 

IT ' ' 



where ^(k, lo) is the spectral function and Qn (k, uj + iS) 
and 033 (k,w -|- i5) are the electron components of the 
superconducting Green's function. Generally, we find: 



A{k,uj) 



ely{2u; -^x- Cy) - (^ + Cy)i^' -e- A?) - (^ + ^x){u^' - ^1 - ^l) 



El -El 

^ -(5(£;3 -u)- S{E3 + uj)) - -^{S{Ei -uj)- 5{Ei + uj)) 



where Ei and E^ are the positive eigenvalues of the matrix -B(k) in ^ (see Ref. [l|). For the case where Ai = A2 = A 
(valid except for the dx2-y2 pairing symmetry), we have the simplified form: 



i?_(k), 



^(k, u;) =r_r-r:i[<5(ijA(k) _ ^) _ s{E^{k) + cu)] + ^i-L|i^[5(ijA(k) _ ^) _ s{E^ik) + c)] (7) 



2E^{k) 



2Ef{k) 



with £'^(k) = ^El{k) + A2(k). This resembles two in- dispersions E±. 
dependent single-band superconductors with the energy 
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FIG. 3: Tunneling differential conductance dl/dV oc 
— / J\f{uj)n'p{Lj — eV) as a function of bias voltage eV mea- 
sured with respect to the Fermi energy, where the temperature 
fcflT = 0.005, the chemical potential /i — 1.6, and the gap size 
Ao = 0.1 for all the different pairing symmetries. 



The spectral function at the Fermi energy A(k, uj = 0) 
contains information about the nodal structure for each 
pairing symmetry, as shown in Fig. [21 The Sx2_i_y2 pair- 
ing symmetry exhibits nodes on the Fermi surface for 
all dopings when = (±7r — ky), (±7r 4- ky) and thus 
only the hole Fermi pockets are fully gapped. The d^y 
pairing symmetry also has nodes for all doping, but in 
this case they occur when k^^y = 0, ±7r and so all of the 
Fermi surfaces are gapless. The d^2_y2 pairing symmetry 
(not shown) exhibits nodes on the Fermi surface of the 
hole pockets for any doping and it has a similar effect 
on the electron pockets as the Sx2y2 pairing symmetry 
which is the dominant pairing symmetry that we found 
m Ref.Q. The Sx2y2 pairing only has nodes on the Fermi 
surface above a critical doping fi 2 since the zeros of 
the gap lie at kx.y = ±7r/2. For /i < 2, the electron 
Fermi surfaces are fully gapped, like the hole Fermi sur- 
faces. In principle, information about the form of the 
Sx2y2 gap can be obtained through ARPES. In the folded 
Brillouin Zone, there are two hole pockets at the F point. 
A cos(fca;) • cos(fcj,) order parameter predicts a larger gap 
for the smaller hole Fermi surface and smaller gap for 
the larger hole Fermi surface. 

Tunneling measurements access the local DOS to a first 
approximation. Specifically, if we assume that both the 
tunneling matrix element and the probe DOS are momen- 
tum independent, then the tunneling differential conduc- 
tance is^i: 



dl_ 

dV 



JV{uj)n'p{uj - eV) 



(8) 



where eV is the bias voltage of the tunneling probe and 
n'p{E) = dnp{E) / dE is the derivative of the Fermi func- 
tion. In the limit of zero temperature, we obviously re- 
cover the DOS. From Fig. [31 we see that a fully-gapped 




FIG. 4: Tunneling differential conductance dl/dV as a func- 
tion of bias voltage eV for the pairing symmetry s^2y2 at 
different dopings. Like Fig. [S] ksT = 0.005 and Ao = 0.1. 



Fermi surface yields a corresponding gap in dl/dV at 
low energies, while for a gapless Fermi surface, the DOS 
grows quasi-linearly with energy at w = 0, a text-book 
result. In particular, Sx2^y2 {d^2_y2) pairing produces 
a four-peak structure in the differential conductance be- 
cause in this case each band sees a different order param- 
eter: the hole (electron) Fermi surfaces are fully gapped 
while the electron (hole) Fermi surfaces are gapless. 

Focussing on Sx2y2 pairing^ (Fig. [¥]), we find that 
the differential conductance smoothly evolves from fully 
gapped to gapless behavior with increasing doping, as 
expected. Moreover, when the doping is large, we ob- 
tain a four-peak structure similar to Sj.2_|_y2 pairing, be- 
cause we also have a fully-gapped hole Fermi surface and 
a partially-gapped electron Fermi surface. While it is 
likely that the material cannot be doped high enough so 
that the Sx2y2 = cos{kx)-cos{ky) superconductor becomes 
gapless (the material will most likely exit the supercon- 
ducting state at such high dopings), we believe that the 
predictions above, in particular the evolution of the dif- 
ferential conductance with doping, could be used in care- 
ful experiments to falsify this order parameter. 

Penetration depth - Measurements of the penetration 
depth in the Fe-based superconductors were the first to 
suggest that the Fermi surfaces are fully gapped2ii2^i2£. 
The experiments show an exponentially temperature de- 
cay of SX{T) = A(T) - A(0). Among the diff'erent order 
parameters studied here, such a scenario is only consis- 
tent with Sx2y2 symmetry at low doping or a constant s- 
wave gap. We now obtain the penetration depth for the 
bare two band superconductor with generic A 1,2 gaps. 

To obtain the penetration depth, we perform a text- 
book exercise. We write the FeAs model in real space 
and introduce a gauge field via the Peierls substitution 
cl^Cj^p — > cl^exp{i j!' A ■ dl)cj,f3 where a, (3 are the two 

orbital indices. We pick a Landau gauge A — Ax and ex- 
pand to second order in A, thus obtaining H{A). The sec- 
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ond order term in A is the diamagnetic current while the 
first order term gives the paramagnetic current, whose 
response must be calculated in linear response. We have: 

HiA) « H{0) - 5](JS(^)A.(^) + ^jf (9) 

i 

hence 

Mi) = -^f^ = m+J'A^)A,{^) (10) 

Using translational invariance, the expectation value of 
the diamagnetic current in the ground state is: 

i 

k J, 

+ ^^(cLck,2 + C^Ck,l) (11) 

where the expectation values of the above operators are 
computed in the appropriate ground state. The param- 



agnetic current is obtained through a correlation function 
in linear response: J^(q,cj) = Q^.^{q,uj)A^{q,u!): 

QtM,^'^n) = ^ j^^dre-"-(jP(q,T)j^(-q,0)) (12) 

This is the vacuum polarization. For the FeAs metal {not 
the superconductor), this is explicitly given by: 

Qxx{(i, ivn) = -TFg X! ^ 

k,m 

Tr (j,(k)G(ic^„ + wn, k + |) Jx(k)G(ic^„„ k - |)) 

(13) 

where ujm — {2m + 1)ttT is a fermionic Matsubara fre- 
quency while Vn — 2m:T is a bosonic one. is the cur- 
rent operator, which is expressed as dH/dk^ in the metal. 
For the response to a magnetic field, the limit that has 
to be taken is, upon analytic continuation, w„ —^Lu + iS, 
= 0, q — + 0. The opposite limit w — > 0, q = gives 
the response to an electric field, and hence the electrical 
conductivity. After tedious but straightforward algebra, 
we obtain for the FeAs metal: 



2 v-^ fdE+Y dn{E+) 
,(q^O,^ = 0) = --E' + ^ ^ +^ 

k 



dE. 



(dEVdn{E_) 8{n{E+) - n{E_)) 



dE^ 



r 



{E+ - E^ 



de_ 



dkrr 



The overall factor of 2 reflects the spin multiplicity. Be- 
sides the usual paramagnetic expression (first two terms 
in Eq. (fT4|) ). the cross-orbital exchange introduces an ex- 
tra second term. We have checked that this paramagnetic 
term completely cancels the diamagnetic ground state ex- 
pectation value, as is required for a metal. We performed 
the same calculation in the superconductor. The charge 
matrix operator in our superconductor is: 



Jn = 



/ 1 

0-100 
10 

Vo 0-1 



(15) 



The current operator uses only the kinetic part of the 



kinetic Hamiltonian and is obtained from the continuity 
equation, giving: 



-^^-2 dkx 



(16) 



where {, } is the anticommutator. The penetration depth 
5\{T) = A(T) — A(0) is proportional to the current- 
current correlation function which uses the Green's func- 
tion of the superconductor, not written here due to space 
restrictions. For the case where Ai = A2, we can write 
the current-current correlation function as: 



Q,,(<7^ 0,LJ = 0) 



k 



fdEj, 



^^xy 

dkx 



dE_ 



{2nF{E'^) - 1 



n'AE^) 



E^ 



^-{2n,{E^)-l)^+'' 



E^ 



(17) 
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FIG. 5: Penetration depth 5\{T) = A(r) - A(0) oc Q:rx{(l 
0, w = 0) close to zero temperature for different pairing sym- 
metries at doping fi = 1.6 and gap size Ao = 0.1. The dxy 
curve has been reduced by a factor of two for clarity. The 
d^2_y2 pairing symmetry (not shown) will have a similar low 
temperature behavior to the d^y and s^2^y2 curves. 



We see that the cross-orbital exchange introduces an ex- 
tra term, similar to the case of the FeAs metal, but 
the largest contribution to the temperature dependence 
arises from the first term. We have obtained the expres- 
sion of the current-current correlation function for gen- 
eral Ai 7^ A2, but we do not include it for space reasons. 

We now plot the low temperature dependence of the 
penetration depth SX{T) — X{T) — A(0) for different su- 
perconducting gaps (see Fig. [5]). As expected, the nodal 
order parameters exhibit a linear T dependence (in the 
absence of impurities) while the nodeless order parame- 
ters exhibit an exponentially decaying penetration depth. 
However, as shown in Fig. [SI one qualitative feature is 
that the cos{kx) ■ cos(fcy) order parameter exhibits, upon 
doping, a weakened exponential decay, a signature that 
the gap on the electron (hole) surface decreases upon elec- 
tron (hole) doping. This is a direct consequence of the 
existence of a line of zeroes in between the electron and 
hole pockets. Above some critical doping, the exponen- 
tial decay of SX{T) in the cos(A:2;) •cos(fcj,) superconductor 
becomes linear (Fig. |6]), a sign that the superconductor 
has become gapless. 

NMR Spin Relaxation Rate and the Coherence Peak - 
Existing experimental results for the NMR spin relax- 
ation time Ti at first sight suggest a d-wave symmetry 
for the order parameter, because there is no coherence 
peak in 1/Ti at Tc and 1/Ti scales like just below 
2^ ^28,29,30,31 ^ These results pose a big challenge for the 
s-wave pairing symmetry or any other nodeless order pa- 
rameter. In the case of a cos(fc2:)'cos(fcy) order parameter, 
although we find that the coherence peak due to inter- 
band contributions is non-existent, the intra-band con- 
tributions still give a coherence peak, although smaller 
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0.2 



FIG. 6: Penetration depth S\{T) = A(r)-A(0) for the pairing 
symmetry s^2y2 at different dopings, where Ao = 0.1. 



and flatter than in a pure s-wave scenario. Neglecting 
the intra-band contributions (which could be justified if 
the broadenings of the inter- and intra-band contribu- 
tions are different) can then explain the observed lack 
of the coherence peak, but in general a small coherence 
peak should be seen in cleaner samples. 

The NMR measurements have been performed on dif- 
ferent atoms in the pnictides, including ^^F and ^^As. 
Experimentally, there is no major difference between the 
1/Ti results on these two atoms. This also poses a chal- 
lenge to the NMR theories because the structure factors 
for F and As are different: while the structure factor for F 
is roughly isotropic in the transferred momentum q, the 
As structure factor is roughly ^(q) — cos{qx/'2) cos{qy/2) 
due to the placement of the As atoms in the center of the 
Fe unit cell (although the As are out of plane, we believe 
the cos((7a;/2) cos{qy/2) faithfully represents the structure 
factor). Hence, for small Fermi electron and hole pock- 
ets, the As NMR measurements should not be sensitive to 
the inter-band contributions, whose transfer wavevector 
(tt, 0) is suppressed by the structure factor. 

The NMR spin relaxation rate at temperature T is 
defined as: 



R 



1 



± lini ^[-^+-(^°)] 
27r ojo^o luq 



where 



/r+_(c.o) =5]A(q)e+-(q,c.o) 



(18) 



(19) 



(q, ujq) is the spin susceptibility in the superconduct- 
ing state and A{q) is the structure factor. Since we are 
dealing with singlet superconductivity we have: 



(20) 
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where is now much simpler due to the fact that the 
spin matrix in a superconductor is the identity matrix: 

a;7i,ki ,k2 

X Tr[G(ki, + c^o))G(k2, zc^„)] (21) 



After Matsubara sums, analytic continuation, and taking 
the imaginary part, for the pure gap case Ai = A2 = A, 
we obtain the following formula for l/{TiT), 



ki,k2 



A(ki)A(k2)+S+(ki)S+(k2) 



-2 1 



A(ki)A(k2)+g-(ki)^_(k2 

i?^(ki)2 
A(ki)A(k2) + £;+(ki)£;-(k2) 

i?^(ki)2 



dn 



dn 



^^.(i.-(k2)-i.-(k,)) 
-6iE^ik,)-E^ik,)) 



(22) 



The first two terms in Eq. (j22p represent the intra- 
band contribution and the third term represents the 
inter-band contribution, which is a contribution between 
the electron and hole pockets. Following Bulut and 
Scalapinor^ we phenomenologically take disorder into 
consideration by broadening the Kronecker delta func- 
tions, e.g. 7rS{E^{k2) - Ef{ki)) = -r/((£;^(k2) - 
i?^(ki)^ + r^). This simple inclusion of disorder works 
well towards explaining the experimental data in the 
cuprate case, and merely serves as a cutoff for the sin- 
gularities in the density of states. We perform the mo- 
mentum integrals by Monte Carlo evaluation: this is nec- 
essary due to the fact that we keep the strong-coupling 
superconductivity and do not make the usual approxi- 
mation which transforms the 4 momentum integrals and 
the delta function into an easy one dimensional integral 
over energies close to the Fermi surface. 

The inter-band and intra-band contributions have dif- 
ferent behaviors as a function of temperature. Owing to 
the fact that for ki on the hole Fermi surface and k2 on 
the electron Fermi surface, A(ki) > while A(k2) < 0, 
we expect the inter-band contribution to lack a coher- 
ence peak around the superconducting transition tem- 
perature, which is indeed what we find below. 

We first consider a uniform structure factor, i.e. 
A{q) — 1. In Fig. [3 we contrast the inter-band con- 
tribution for the 3^2 y2 pairing symmetry with that of its 
absolute value, i.e. | cos(fc2:) •cos(fcj^)|, which does not ex- 
hibit a sign change between the hole and electron pock- 
ets. Clearly, the former case does not possess a coherence 
peak, while the latter does, as expected. In Fig. [U we 
plot the intra-band contribution and the total l/Ti for 
both cases. We see that, compared to the absolute value 
case, the coherence peak in 1/Ti is suppressed in the 
cos(A;^) • cos(fcy) case. 

Using the structure factor A(q) for As atoms (Fig. [5]) 
we find that the inter-band component of the to- 




FIG. 7: Monte-Carlo calculation of the (normalized) 
inter-band contributions to the NMR coherence peak for 
Aocos(fca;) ■ cos(fcj;) (red) and a fixed-sign version of it, 
Aolcos(fc2,) • cos(fcy)| (blue). We choose a large Ao — \ti\/5 
and Ao/Tc = 2. The broadening factor is F = Tc/S, and 
fi = 1.8 corresponding to 18% electron doping. Inset: Tem- 
perature dependence of the NMR spin relaxation time for 
Ao cos(fca;) cos(fcy) (red). The structure factor here is taken 
to be A{q) = 1. 



tal NMR spin relaxation rate decreases. While for 
A(q) = 1 the inter-band contribution represented about 
1/6 of the overall spin relaxation rate, for A{q) — 
cos((7a;/2) cos((7j,/2) that ratio decreases to about 1/12. 
We hence find that the intra-band contribution is domi- 
nant in the case of the As structure factor. However, we 
also find that the structure factor reduces the intra-band 
coherence peak, to give an overall result plotted in Fig. [9l 
Finally, we find that the NMR relaxation rates for 
the nodal superconductors dxy and s^^+yi, depicted in 
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FIG. 8: Monte-Carlo calculation of the (normalized) intra- 
band contributions to the NMR coherence peak for the 
Aocos(fca;) • cos(fcj;) (green) gap (the intra-band contribution 
is equal for the two gaps Aq cos(fca;) •cos(fcy) and Ao| cos(fca;) ■ 
cos(fcj;)|). The total, intra plus inter band contributions for 
Ao cos{kx) ■ cos{ky) (red) and Ao| cos{kx) ■ cos(fcj,)| (blue) are 
also plotted. We can see that the intra-band contribution is 
hence much larger than the inter-band contribution for both 
these order parameters, and hence the Aocos(fci) ■ cos(/cy) 
gap should exhibit a small coherence peak. We choose a 
large Aq — \ti\/5 and Aq/Tc = 2. The broadening factor 
is r = Tc/5, and /i — 1.8 corresponding to 18% electron dop- 
ing. The structure factor here is taken to he A{q) = 1. 



FIG. 9: Monte-Carlo calculation of the (normalized) intra- 
band contributions to the NMR coherence peak for the 
Aocos(A:j;) • cos{ky) (green) gap, with the structure factor 
yl(q) — cos{qx/2) cos(gy/2) for the As atoms. The total, intra 
plus inter band contributions for Aocos(fc2,) • cos(fcy) (blue) 
are also plotted. The coherence peak is diminished from the 
case when A(q) — 1, plotted previously. Inset: Tempera- 
ture dependence of the inter band contribution. We choose a 
large Aq = |ii|/5 and Aq/Tc ~ 2. The broadening factor is 
F = Tc/5, and fi = 1.8 corresponding to 18% electron doping. 



Fig. [ini lack a coherence peak as expected. 

We predict that future experiments will see a small 
coherence peak resulting from the intra-band contribu- 
tion. Our results show that, barring different scatter- 
ing rates for inter- and intra-band scattering, the overall 
intra-band contribution to the NMR relaxation rate is 
roughly a factor of 5 times larger than the inter-band con- 
tribution. This can also be argued on general grounds, 
provided that the hypothesis of weak-coupling theories 
and LDA (i.e. there is a quasi-nesting of the electron 
and hole Fermi surfaces in the parent material) is cor- 
rect. Upon doping with either electrons or holes, either 
the electron or hole Fermi surfaces will become consid- 
erably larger than the other one. This means that the 
inter-band contribution to the NMR spin relaxation rate 
diminishes: it of course vanishes if one could, theoreti- 
cally, deplete one of the Fermi pockets. Meanwhile, the 
intra-band contribution should, on general grounds, re- 
main roughly constant upon doping because the overall 
size of the sum of the Fermi surfaces is relatively con- 
stant. All these general arguments are supported by our 
explicit calculation. 

A few other remarks about the NMR spin relaxation 
rates are in order: (i) The observed temperature de- 
pendence of 1/Ti cannot be viewed as evidence against 
s-wave pairing symmetries. In fact, the temperature 
dependence just below Tc is very sensitive to the ratio 



0.6 O.S 1.0 



0.2 0.4 0.6 0.8 



FIG. 10: Monte-Carlo calculation of the (normalized) intra- 
band (green) and total (blue) contributions to the NMR co- 
herence peak for the Ao sm{kx) sin(fci,) (left) and Ao(cos(fca,) + 
cos(fcj;)) gap. These are nodal superconductors and lack a co- 
herence peak. We choose a large Aq — \ti\/5 and Aq/Tc ~ 2. 
The broadening factor is F = Tc/5, and fi — 1.8 correspond- 
ing to 18% electron doping. The structure factor here is taken 
to be A{q) = 1. 



A/ksTc. We find that the behavior can be obtained 
by choosing A/ksTc ~ 2 for our large gap value, and the 
power of the temperature dependence can increase even 
further by increasing this ratio, (ii) Although we predict 
that there should be a coherence peak in the clean limit, 
impurities can efficiently reduce the coherence peak in 
a two-band system. A weak inter-band impurity scat- 
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tering but strong intra-band scattering can suppress the 
coherence peak. This has been investigated in MgB2^ 
where the coherence peak is also not easily observed ex- 
perimentally"^*'. Since the superconductivity in Fe-based 
superconductors is created by doping, it is reasonable to 
assume that disorder is stronger than that in MgB2. To 
observe the coherence peak, we require a very clean sam- 
ple, (iii) Our calculation is based on a two-band model. 
This model can be over-simplified when one tries to use it 
to predict quantitative experimental measurements. For 
example, the detailed shape of Fermi surfaces and its 
doping dependence may not be quantitatively accurate. 
Therefore, the predictions in this paper with regard to 
doping concentration should be viewed as qualitative. 

Conclusion - We have calculated the spectral func- 
tion, the DOS, the tunneling differential conductance, the 
penetration depth, and the NMR spin relaxation rate for 
different superconducting order parameters in the iron- 
pnictides. We have emphasized that the nodal structure 
of the Sx2y2 order parameter will result in a qualitative 
change in these experimental observables with increasing 
doping, as the superconductor crosses over from gapped 
to gapless. Thus, one can in principle probe the exis- 
tence of this pairing symmetry in the iron-pnictides by 
analyzing the behavior of the spectral function, the DOS 



and the penetration depth as a function of doping. For 
the 1/Ti NMR spin relaxation rate, if only the inter- 
band contribution is considered, our theoretical results 
are consistent with the current experimental results, in- 
cluding the T"^ behavior and the absence of a coherence 
peak. However, by including the intra-band contribution, 
a small coherence peak at the transition temperature will 
be present in a clean sample although it is smaller than 
that in a sign-unchanged s-wave. 

Note - During the completion of this work, we be- 
came aware of two recent papers that also calculate the 
spin-lattice relaxation rate for the 5^,2^2 order parameter 
in the iron-pnictide o— , and another recent paper that 
considers the experimental consequences of two different 
pairing symmetries^!. 
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